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ABSTRACT 

We examine the conditions of the plasma along a sample of 'classical' Herbig-Haro (HH) jets located in the Orion and Vela star forming 
regions, through combined optical-infrared spectral diagnostics. Our sample includes HH 111, HH 34, HH 83, HH 73, HH 24 C/E, HH 24 J, 
observed quasi-simultaneously and in the same manner at moderate spatial/spectral resolution. Once inter-calibrated, the obtained spectra cover 
a wide wavelength range from 0.6 - 2.5 /jm, including many transitions from regions of different excitation conditions. This allows us to probe 
the density and temperature stratification which characterises the cooling zones behind the shock fronts along the jet. From the line ratios we 
derive the variation of the visual extinction along the flow, the electron density and temperature (n e and T e ), the hydrogen ionisation fraction x e , 
and the total density n H in the emission region of different lines. The knowledge of such parameters is essential for testing existing jet models 
and for planning follow-up high-angular resolution observations. 

From the diagnostics of optical forbidden lines we find, on average, that in the examined jets, in the region of optical emission, n e varies 
between 50 crrT 3 and 3 10 3 cur 3 , x e ranges between 0.03 and 0.6, and the electron temperature T e is ~1.3 10 4 K in the HH 111 and HH 34 jets, 
while it appears to be higher (1.8 10 4 K on average) in the other examined jets. The electron density and temperature derived from [Fen] lines, 
turn out to be, respectively, higher and lower in comparison to those determined from optical lines, in agreement with the fact that the [Fe n] 
lines arise in the more compressed gas located further from the shock front. An even denser component in the jets, with values of n e up to 10 6 
cm -3 is detected using the ratio of Calcium lines. 

The derived physical parameters are used to estimate the depletion onto dust grains of Calcium and Iron with respect to solar abundances. This 
turns out to be quite substantial, being between 70% and 0% for Ca and ~90% for Fe. This leads us to suggest that the weak shocks present in 
the beams are not capable of completely destroying the ambient dust grains, confirming previous theoretical studies. We then derive the mass 
flux rates, Mj ct , in the flows using two independent methods. Taking into account the filling factor of the emitting gas, Mj et is on average 5 
10~ 8 M yr _1 . The associated linear momentum fluxes (Pj et = v ; «Mj ct ) are higher than, or of the same order as, those measured in the coaxial 
molecular flows, where present, suggesting that the flows are jet driven. 

Finally, we discuss differences between jets in our sample. In general, we find that higher ionisation and electron temperatures are associated 
with less dense jets. The comparison suggests that the shock mechanism exciting the knots along the flows has the same efficiency in all the 
examined objects, and the observed differences are consistent with the different densities, and hence cooling rates, found in the various flows. 
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citation properties are reasonably well known. This in principle 
allows us to go back from the observed line intensities and line 
ratios to the physical conditions of the emitting material, which 
in turn provide essential information to test theoretical mod- 
els and plan challenging observational programs, e.g. those in- 
volving interferometry. For quite a long time line diagnostics 
have been used mainly to determine the gas electron density 
(n e ) and electron temperature (T e ) (see, e.g. Osterbrock 1994). 
The spectra of line emitting regions, however, contain much 
more information and the development of a number of new di- 



1. Introduction 

Spectral analysis is a powerful tool to investigate the nature and 
properties of astrophysical nebulae excited by shocks and/or 
energetic radiation. The diffuse matter in such objects produces 
a wealth of permitted and forbidden emission lines, whose ex- 
Send offprint requests to: Linda Podio, lindapod@arcetri.astro.it 
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agnostic techniques allow us to put, directly or indirectly, new 
important constraints on the gas physics. The study of stellar 
jets is one of the fields that has benefitted enormously from 
progress in spectral techniques. These collimated flows appear 
right from the birth of a star and indeed are believed to play 
a central role in the star formation process itself (see, e.g., 
lEisloffel et al. 20001 |ReTpurth & Bally 200 1| |Ray et al 2003] 
IBacciotti et al. 20041 . In fact, jets can help clear the circum- 
stellar environment, thus setting a limit to the final cen- 
tral mass, and may even extract most of the excess angu- 
lar momentum from the star/disk system ( Bacci otti et al. 2 002 
ICoffey et al. 2004| IWoitas et al. 2005> . Prototypical cases are 
the HH 34 and HH 111 jets in Orion ( |Reipurth & Bally 2001) , 
also studied in this work. All the observed features in stel- 
lar jets, such as bright knots along the flow and giant bow 
shocks, are believed to trace shocks that develop within these 
highly supersonic flows. Passing through the shock front, part 
of the bulk kinetic energy of the flow is turned into thermal mo- 
tions, and the gas is suddenly heated to high temperatures (10 5 
K or more), compressed, and partially ionised. This situation 
favours the collisional excitation of upper levels of atomic tran- 
sitions. A number of radiative lines, whose type and strength 
are characteristic of shock excitation, are produced. These, 
in turn, radiatively cool down the gas to its original temper- 
ature. The whole process takes place in a limited region of 
space behind the shock front called the 'cooling zone', whose 
length may vary from 10 13 to about 10 15 cm, depending on 
the shock strength, the pre-shock conditions and the compo- 
sition of the gas. The structure of radiatively cooling shocks 
has been widely studied in the past by several authors (see e.g. 
Hartigan et al. 1987| iHollenbach 19891 |Hartigan et al. 1994 
Hollenbach 1997IIFlower et al. 20031 |Hartigan 2003| and refer- 
ences therein). Within the cooling region individual physical 
quantities vary enormously, creating a stratification of excita- 
tion conditions that results in the production of different kinds 
of lines. In ground based observations at moderate spatial reso- 
lution, such as those presented here, the internal structure of the 
shock cooling zone remains spatially unresolved. Nevertheless, 
the stratification of such regions can be investigated spectrally 
by analysing the many different lines arising from different lay- 
ers. 

In the optical and Near Infrared (NIR) ranges used in this 
work one observes forbidden and permitted lines of abundant 
atomic and molecular species, such as H, He, O, S, N, Fe, C, 
Ca, H2 (Reip urth & Bally 200 l| l. To extract more information 
from these lines, several attempts have been made to use novel 
techniques by various groups. Bohm et al. (1980), Brugel et 
al. (1981) were among the first to examine the potential of 
combining the observed line ratios using low spatial resolu- 
tion spectra integrated over the shock cooling zone. Hirth et 
al. (1997) examined long-slit spectra to study the spatial and 
kinematical properties of the forbidden emission line regions 
and micro-jets of a large sample of T Tauri stars. Intensities of 
atomic lines have been compared with the prediction of shock 
models in, e.g., Raga & Bohm (1986), Hartigan et al. (1987), 
and Hartigan et al. (1994), indirectly inferring quantities such 
as the hydrogen ionisation fraction x e and the intensity of the 
magnetic field in the pre-shock region. A new simple technique 



to measure x e and T e from the ratio of optical lines was first pre- 
sented by Bacciotti, Chiuderi & Oliva (1995), and subsequently 
refined in Bacciotti & Eisloffel (1999) (hereafter BE99). The 
method, referred to as the 'BE' technique, is based on the fact 
that the gas emitting forbidden lines is collisionally excited, 
but no assumption is made regarding the heating agent. This is 
clearly an advantage if the results are used to validate a given 
thermal model. On the other hand, the method assumes that the 
emitting gas is at a single temperature, which is not true in the 
cooling region behind a shock front. In this regard, one should 
consider that the results of the BE technique are relevant to 
the region behind the shocks in which the employed lines have 
their peak emission (see the discussion and diagrams in BE99). 
As we show in this paper, the stratification of temperature and 
densities present in the entire cooling zone can be traced us- 
ing a larger sample of lines appropriate for different excitation 
conditions. 

A number of jets have been analysed with the BE tech- 
nique, leading to the finding that jets are only partially ionised, 
with average x e values between 0.01 and 0.6. This procedure is 
much easier to apply than a grid of shock models (although 
it is of more limited application), and thus well suited for 
the analysis of big datasets, as those provided by high angu- 
lar resolution observations (Ba cciotti, Eisloffel & Ray, 1999| 
Bacciotti 2002 1. The main advantage in determining x e , inde- 
pendently from the gas heating mechanism, is that by combin- 
ing its value with the derived electron density one can estimate 
the total density, n#, a fundamental parameter of the jet that is 
critical in the various models. 

Diagnostic techniques using both optical and NIR lines, on 
the other hand, provides important complementary informa- 
tion on jet parameters. In particular, the combination of op- 
tical and NIR [Fen] lines gives an independent tool to de- 
termine T e and n e , which does not rely on the choice of ele- 
mental abundances, in the more compressed and cooler post- 
shock layer where these lines are excited ( Ni sini et al. 2 002 
IPesenti et al. 20031 |Hartigan et al. 2004| i. In addition, the NIR 
H2 lines provide a means to probe the molecular component 
of shock excited gas, which may give a significant contribu- 
tion to the gas cooling in low velocity, magnetized shocks 
(Eisloffel, Smith and Davi s 2000|lGiannini et al. 2004> . 

Very recently, we have investigated the advantages of a 
combined optical and NIR spectral analysis of the HH 1 jet 
( Nisi ni et al. 20 05 ) using a variety of diagnostic lines over a 
wide wavelength range (from 0.6 to 2.2 fim). The adopted pro- 
cedure turned out to be extremely powerful for constructing 
detailed physical maps of the stratified media in the beams of 
stellar jets. In this paper, we analyse a large sample of 'clas- 
sical' stellar jets, located in the Orion and Vela star formation 
regions, with the same technique. We obtain their basic param- 
eters, as well as important quantities derived from them, such 
as the mass flux rates. We analyse the variations of the parame- 
ters both behind each shock, i.e. depending on the tracers, and 
along the jet, i.e. as a function of the distance from the emit- 
ting source. We then discuss the differences/similarities among 
jets having different properties and ages. This information can 
in turn be used for the selection of suitable candidates for high 
angular resolution observations. 



L. Podio et al.: Stellar jets analysed with opt/NIR diagnostics 



3 



This paper is organized as follows: we describe our obser- 
vations in Sect. 2 and briefly recall our diagnostic procedure in 
Sect. 3, adding a detailed discussion on the choice of the best 
set of elemental abundances, an issue that was not examined in 
our previous papers. In Sect. 4 we describe the results obtained 
for each target in the sample, and we discuss them in Sect. 5 
together with the derivation of the depletion of refractory el- 
ements and the mass and linear momentum flux rates. Sect. 6 
summarizes our findings. 

2. Observations 

We observed a sample of classical protostellar jets (HH 34, 
HH 111, HH 83, HH 73, HH 24 C/E, HH 24 J) in the spec- 
tral range from 6015 A to 2.52 yum. We acquired the optical 
spectra (6015-10320 A) at the ESO 3.6-m telescope equipped 
with the spectrograph EFOSC2. Infrared spectra were taken at 
the NTT with Sofl using both the blue grism (IR-GB: 0.95- 
1.64 fim) and the red grism (IR-GR: 1.53-2.52 //m). In both 
cases the spectra were taken aligning the slit parallel to the jet 
axis. The position angles are: 166° for HH 34, 277° for HH 
111, 299° for HH 83, 303° for HH 73, 334° for HH 24 C/E, 
and 309° for HH 24 J. The observational settings were chosen 
to obtain, as far as possible, homogeneous spectra in the two 
wavelength ranges, in order to be able to apply a combined op- 
tical/NIR analysis. The observations were all carried out over a 
short period of time (7-8 and 11-12 January 2003, for the opti- 
cal and IR observations, respectively). In this reagard, note that 
jet parameters change over timescales not shorter than a few 
years and probably decades in most cases. A 1" width slit was 
used for both the optical and the IR spectra, in order to cover 
the same section of the jet and to have similar spectral reso- 
lution (R ~ 600). Moreover, the spatial scale of the two cam- 
eras is comparable (O'.'314/pixel for EFOSC2 and 07288/pixel 
for Sofl). The integration times for the EFOSC2 observations 
were 1800s for HH 111, HH 34, HH 73 and 5400s for HH 
83, HH 24 C/E, HH 24 J. In the infrared the exposure times 
for each of the two Sofl grisms are: 1800s for HH 111 and 
HH 34, 1200s for HH 83 and HH 24 C/E, 600s for HH 24 J, 
800s in the IR-GB and 400s in the IR-GR for HH 73. In addi- 
tion, telluric and spectro-photometric standards were observed 
in order to correct for the atmospheric spectral response and to 
flux-calibrate the spectra respectively. The wavelength calibra- 
tions were performed using a helium-argon lamp in the optical 
and a xenon lamp in the infrared. The data reduction was done 
using standard IRAF tasks. We obtained for each jet three in- 
dividually calibrated spectra in the optical, IR-GB and IR-GR 
ranges. For each knot, a single spectrum from 6015 A to 2.52 
yum was then formed as follows: (i) The same spatial reference 
was defined in each of the three spectra, (ii) A knot length was 
set by comparing the spatial profiles along different lines with 
our images acquired through [SII] and H2 narrow-band filters 
and with high angular resolution images from the literature. 
The same knot size was then adopted in all three spectral do- 
mains, (iii) We then extracted from each spectral image, and for 
each knot, the corresponding spectrum integrating over the size 
defined above, (iv) For each knot, the three obtained spectra 
were inter-calibrated measuring the fluxes of the lines located 



in the overlapping spectral regions (the [CI]9850A line flux in 
the optical/IR-GB and the [Fell] 1 .64/jm line in the IR-GB/IR- 
GR). 

3. Description of the combined optical/NIR 
spectral analysis 

In Sect. 3.1 the principles underlying our diagnostic techniques 
are briefly recalled. In Sect. 3.2 we present a critical discussion 
on the adoption of a given set of elemental abundances, which 
is necessary for our analysis. 

3.1. Derivation of the physical conditions along the jet 
beam 

Thanks to the many lines observed over the wide spectral range 
covered by our EFOSC2/SofI spectra, and using together opti- 
cal JBE99> and infrared (INisini et al. 2002IIPesenti et al. 2 003 1 
diagnostic techniques, we have been able to determine the pa- 
rameters describing the stratified medium behind each shock 
even if the knots along the jets are not spatially resolved. 

Firstly, the visual extinction, Ay, has to be determined 
to deredden the spectra and be able to combine lines 
that are far apart in wavelength. Ay can be found using 
pairs of distant lines coming from the same upper level 
of a single ion and assuming a reddening law. Suitable 
choices in our case are the [Feii]1.64yKm/1.25/im or the 
[Fen]1.64//m/1.32yum ratios. We tried to calculate Ay from 
both these ratios and using different sets of Einstein co- 
efficients (Nussbaumer & Storey 1988, IQuinetetal. 1996| 
|Smith & Hartigan"2 006l. With any of the above sets the 
[Fen]1.64 / um/1.25 / um ratios give higher values of Ay than 
those inferred through using [Fen]1.64jum/1.32yum. Nisini et 
al. (2005) pointed out that the [Fen]1.64/im/1.25/im ratio, 
compared with the theoretical value adopting the radiative rates 
from Nussbaumer & Storey (1988) provides extinction values 
which are too high, inconsistent with determinations in the op- 
tical and with model predictions (see their Appendix B). We 
find that the [Fen]1.64yum/1.25/um ratio empirically estimated 
by Smith & Hartigan (2006) from the spectrum of P Cygni 
produces even higher Ay values. On the other hand, Smith 
& Hartigan (2006) derive a [Fen]1.64 / um/1.32 / um theoreti- 
cal ratio very similar to that of |Nussbaumer & Storey (1988) 
which gives an Ay more consistent with other independent 
determinations (Nisini et al. 2005). Thus we adopted the 
Ay value derived from the [Fen]1.64yum/1.32yL(m ratio and 
the [Nussbau mer & Storey (1988)] coefficients as in Nisini et 
al. (2005). All the spectra were then corrected for extinc- 
tion using a standard dereddening procedure dDraine 1 989 1 
and an interpolation of the extinction law derived by 
|Rieke & Lebofsky (1985)| for the near-IR bands. 

The electron density, n e , in the region of optical emis- 
sion, was derived from the [S n]/l6716A/[S ii]/16731A ra- 
tio (Osterbrock 1994). Since in some of the knots the 
[S n]/l/l6716, 6731 doublet was not resolved sufficiently to 
separately measure the line fluxes, we used data taken two 
years before by our group using EFOSC2 at higher spec- 
tral resolution (Medves, Bacciotti & Eisloffel, in prep.). 
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Intervening proper motions of the jet knots can be neglected 
for the purposes of our analysis. Then, using the BE tech- 
nique (see BE99), that employs selected optical transitions 
of S + , O , and N + , we derived the electron temperature, 
T e , and the hydrogen ionisation fraction, x e . Since in de- 
termining these quantities one uses ratios between different 
species (namely [N n]^(6548+6583)/[Oi]^(6300+6363) and 
[Sn]Xl(6716+6731)/[Oi],Ll(6300+6363)), the technique re- 
quires the adoption of a given set of elemental abundances. 
This issue is examined in detail in the next Section. Finally, 
from the values inferred for n e and x e , we derived a gross es- 
timate of the total hydrogen density (nH-n e Jx e ). Refinements 
to this density estimate are discussed in Sect. 5. Also, as men- 
tioned in the Introduction, we stress that the results of the BE 
technique are only relevant to the region of the cooling zone 
behind the shock front where the considered optical lines emit 
the most (see BE99l. 

The errors that affect the parameters obtained through the 
BE diagnostic are due to measurement errors of the line fluxes 
(which depend on the signal-to-noise ratio) and the uncertainty 
in the determination of Ay which is used to deredden the line 
fluxes. To determine x e and T e the uncertainty in n e values also 
has to be taken into account. On the other hand, since the lines 
used in this diagnostic are very close in wavelength the main 
source of uncertainty is the measurement error. In general, for 
the brightest knots, where the signal-to-noise is high (S/N > 
10), the errors are <5% for n e , <15% for x e values and <10% 
for T e . 

[Fen] lines give a completely independent way to de- 
termine the physical conditions of the gas. From the ratios 
[Fen]1.64/vm/1.53yum and [Fen]1.64/mT/1.60/mi one can in- 
fer the electron density n e in the region of Fe emission. The 
ratios between optical and infrared [Fen] lines can be used 
to find the temperature in the same region. The relevant ra- 
tios for this goal, are those of the 1.64/mi line with different 
transitions between 8000 A and ljt/m originating from the a 4 P 
term dNisini et al. 2 002 Pes enti et al. 2 003). As we will show, 
[Fe n] lines trace the cooler and denser gas located further from 
the shock front than the zone of optical line emission, and thus 
probe a region where the emitting material is more compressed 
(see lHollenbach 1997i . 

In order to check if there is an even denser compo- 
nent in the knots we used the Caii/l8540/[Caii]/l7290 and 
[Fe ii]/17155//18617 ratios. The theoretical values of Ca n ratios 
for three different values of the temperature (T = 5000 K, T = 
10000 K, T = 15000 K) were computed with a five-level sta- 
tistical equilibrium code (radiative transition rates from NIST; 
collisional rates from Mendoza 1983 and Chidichimo 1981 1. 
This code assumes purely collisional excitation. The observed 
ratio [Can]/l7290//i7324, ~1.5 in all the cases, is consistent 
with the ratio expected for collisional excitation. Such a ratio 
is however also consistent with fluorescence pumping of the 4p 
levels giving rise to the Ca n H and K UV lines, followed by 
cascade to the 3d level ( Hartig an et al. 2004| l. Hartigan et al. 
(2004) estimated a neglegible fluorescence pumping rate with 
respect to the collisional rate in the Ca n excitation of the HN 
Tau jet observed at 10 AU from the source. We assume that this 
is the case also for the jets in our sample, since we are deal- 



Table 1. Abundance sets 



Abundance Set 


S/H 


O/H 


N/H 




(io- 5 ) 


(10- 4 ) 


(io- 5 ) 


(i) Solar (G&S 98)" 


2.14 


6.76 


8.32 


(ii) Solar (A 05)* 


1.38 


4.57 


6.03 


(iii) Adopted in BE99 1 


1.6 


6.0 


11.0 


(iv) Orion (E 04) rf 


1.66 


4.47 


5.37 







IGrevesse & Sauval 1 998 
* |Asplund et al. 2005| 
c Bacciotti & Eisloffel 1999 
lEsteban et al. 20041 



ing with jets located at large distance from the driving source, 
where no significant UV field is expected. The observed val- 
ues of [Fen] ratios, instead, were compared with prediction 
diagrams of Bautista & Pradhan (1996) in order to determine 
the electron density. These diagrams plot the predicted [Fe n] 
ratios versus the electron density for a fixed temperature (we 
considered the curve for T = 10 4 K). 

Finally, to derive the gas physical conditions in the region 
of molecular emission, we use H2 lines diagnostics. In par- 
ticular, from the rotational diagrams we estimate the temper- 
ature and the total column density of the molecular gas (see 
INisini et al. 200 2 for more details). 

3.2. Choice of elemental abundances 

As discussed in Sect. 3.1, to estimate the ionisation fraction, 
x e , and the temperature, T e , one has to assume the relative 
abundances of S, O, N used in the diagnostic. As a conse- 
quence, the derived values may depend on the adopted set of 
abundances, an issue that has to be investigated in more de- 
tail. For example, in Nisini et al. (2005) we performed our di- 
agnostic analysis adopting solar abundances from Grevesse & 
Sauval (1998). The majority of our targets, however, are lo- 
cated in the Orion Cloud. Moreover, there are many lines of 
evidence that suggest Solar System abundances may not even 
be representative of the local ISM 4.6 Gyr ago, at the time the 
Sun formed (Wilson & Roo d 19941 . It follows that solar abun- 
dances may not be adequate to study the properties of the jets 
in our sample. To elucidate the situation, we test the effect of 
abundance variation on our diagnostic results. The abundance 
sets we consider are: (i) the solar abundances determined by 
IGrevesse & Sauval 19981 (G&S 98), (ii) the most recent de- 
terminations of the solar abundances by Aspl und et al. 2005] 
(A 05), (iii) the abundances adopted in BE99] (taken from 
|Hartigan et al. 1987) with the aim of testing agreement be- 
tween the BE diagnostics and shock models and (iv) the abun- 
dances determined for the interstellar gas in the Orion Nebula 
by Esteb an et al. 20 04 (E 04). The values of elemental abun- 
dances of S, O, N in these sets are summarized in Tab.^ 

We computed the values of x c and T e for the HH 111 jet 
assuming the above sets of abundances. The results are shown 
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HH 111: ABUNDANCE TESTS 




0.1 



30 50 70 

distance from source (orcsec) 



A-B E F G21 1 H II J K1 L M N 



T-y 



♦ observed ratios 

G&S 98 ~ 

- - - BE 99 : 
E 04 

A 05 





10 30 50 70 145 155 

distance from source (arcsec) 

Fig. 1. Variations in results using the BE technique and 
different adopted elemental abundance sets. The HH 1 1 1 jet is 
used for test purposes. Upper panel. Values of the ionisation 
fraction x e and electron temperature T e in each jet knot 
assuming: solar abundances from IGrevesse & Sauva~1 998 
(squares) and Asplu ncTet al. 2005] (crosses); abundances 
adopted in IBE99I (triangles'): Orion abundances from 
Este ban et al. 20041 (circles). Lower panel: Comparison be- 
tween observed values of the [Ni](1.04;um)/[Nii](6583A) and 
[Sn](1.02+1.03jum)/[Sn](6716+6731A) ratios and theoretical 
predictions for the various abundance sets. 

in the upper panels of Fig. ^ The figure shows that with the 
Orion abundances of Esteban et al. (2004) we obtain the high- 
est values of x e and T e from our analysis. On the other hand, 
these values are within 15% of those obtained using the recent 
solar abundances determination of Asplund et al. (2005), since 



the abundances of the elements involved in the diagnostics are 
quite similar in the two sets. Moreover, the values of x e and 
T e found assuming the solar abundances by Grevesse & Sauval 
(1998) are again in good agreement with the previous ones. 
The reason here is somewhat different: despite the Orion E 04 
set and the solar G&S 98 set having different abundances, the 
ratios O/N and S/O used in the analysis turn out to be similar. 
On the other hand, the values inferred from the abundance set 
adopted in BE99 show differences of up to 40% for x e and up to 
25% for T e . We also tested the results obtained with the differ- 
ent sets using the observed ratios [Ni](1.04^m)/[Nn](6583A) 
and [Sii](1.02+1.03jum)/[Sn](6716+6731A), that do not de- 
pend on abundances, but do depend on the derived x e and T c . 
Comparing observed and theoretical values (see the lower pan- 
els of Fig. ^ we find that the abundance sets (i), (ii) and (iv) 
are all consistent with the observations, while the abundances 
adopted in BE99 that yield the lowest theoretical values for 
these ratios, deviate more. 

Following the above results, we chose to adopt the abun- 
dance set by Esteban et al. (2004) for the jets in our sample 
located in the Orion cloud. For HH 73, located in Vela, we ten- 
tatively assume the same, since no abundance determination 
exists for Vela. 



4. Diagnostic results: basic physical parameters of 
the examined jets 

In this Section we describe our determination of basic physical 
parameters, i.e., density, temperature and ionisation fraction, in 
each jet knot, obtained from the application of the combined 
optical/NIR technique to our sample. The best data were ob- 
tained for the HH 111 jet, allowing us to get very accurate val- 
ues in this case. Parameters for HH 34 also appear to be very 
reliable. HH 111 and HH 34 are very well known HH objects 
and they have been investigated previously both in the optical 
wavelength range (Eisloffel & Mundt 1992, Morse et al. 1993, 
Noriega-Crespo et al. 1993, Hartigan et al. 1994, Eisloffel & 
Mundt 1997, Reipurth et al. 1997, |BE^9l Hartigan et al. 2001, 
Reipurth et al. 2002, Raga et al. 2002) and in the near-infrared 
(Stapelfeldt et al. 1991, Stanke et al. 1998, Davis et al. 2001, 
Nisini et al. 2002). Our purpose here is to combine the infor- 
mation coming from the different tracers emitting in the two 
wavelength ranges, in order to obtain a detailed physical map 
of the stratified medium in each unresolved cooling zone along 
these outflows. The fainter objects HH 83, HH 73, and HH 24, 
instead, were either not seen in the IR (HH 73) or were visible 
only in one or two knots (HH 83 and HH 24). Thus for these 
jets we only apply the optical diagnostic to obtaine a descrip- 
tion of the physical conditions of the gas emitting in this range. 
Note that, although it is not possible to estimate the extinction 
when the [Fe n] lines are not detected, we can still use the opti- 
cal diagnostic since the lines involved in the BE technique are 
very near in wavelength. A summary of our results is given in 
Tab. |2 but note that the derived physical parameters are aver- 
ages for the brightest knots. 
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Fig. 2. Spatial profiles of the [Sn], [Fen], and H2 2.12 fim 
lines along the HH 111 jet. The zero point of the spatial scale 
is the source HH 111 IRS (RA(2000): 05 51 46.1, Dec(2000): 
+02 48 30). 



4.1. HH111 

One of the best known stellar jets is HH 111. This outflow, lo- 
cated in the L1617 cloud in Orion (D = 460 pc), is powered by 
the young star IRAS 0549 1 +0247 and extends over several par- 
sees ( |Reipurth et al. 1997) . The exciting source is deeply em- 
bedded in the parental molecular cloud core (in our spectra the 
source is detected only in the infrared K-band) and the base 
of the jet is not visible at optical wavelengths. The main visi- 
ble jet emerges from the cloud only 15" from the source and 
is blue-shifted. Adopting the nomenclature of Reipurth et al. 
(1997), this lobe consists of a long chain of knots, with a bow- 
like morphology, observable in the optical from 15" (knot A-B) 
to 80" (knot P) and terminates with a bow shock (knot T-V) at a 
distance of 150" from the driving source. The red-shifted lobe 
is almost completely obscured, with only knot ZL visible in the 
optical, while knots ZO and ZV (symmetric to knots O and T-V 
in the blue lobe) have been observed only in the near-infrared. 

Several authors have investigated the physical properties of 
this flow in the past. Morse et al. (1993) measured the elec- 
tron density along the jet from the [SII] line ratio and presented 
a complete list of the lines emitted in the optical range (see 
also |Noriega-Crespo et al. 1993| l. Hartigan et al. (1994) esti- 
mated a value of the ionisation fraction averaged along the en- 
tire flow comparing the observed line ratios with the predic- 
tion of shock models, and derived an estimate of the average 
mass loss rate. Internal gas kinematics has been investigated 
by Hartigan et al. (2001) and Raga et al. (2002). In particular, 
in the first paper proper motions of the knots were determined, 
from which shock velocities of about 40 km s were derived, in 
agreement with the type of line observed. In the infrared wave- 
length range, Davis et al. (2001) studied the kinematics of the 
gas emitting in H2, while an analysis of H2 and [Fen] lines 
along the jet has been presented in Nisini et al. (2002). 

In this paper we apply our combined optical/NIR analysis 
to the knots in the blue lobe from 15" outwards, where both 
optical and infrared lines are observed (except for knot T-V 
that is out of the frame in the IR spectra). In Fig. |2 the spatial 
normalized distribution of the optical [S n] line is compared 
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Fig. 3. Variation of the derived physical parameters along the 
HH 111 jet. From top to bottom panel: intensity profiles of the 
optical lines, the electron density, n e , in units of 10 3 cm -3 , the 
ionisation fraction, x e , the temperature, T f , in units of 10 4 K 
and the total density, n#, in units of 10 4 crrT 3 . The open circles 
are the values derived from the [Fen] lines, while the filled 
circles are parameters inferred from the optical S + ,N + and 0° 
lines using the BE technique. Note that the zero point of the 
spatial scale is the driving source HH 111 IRS (see Fig|2J. 

with the corresponding [Fen] and H2 2.12 jjm distributions. 
[S n] and [Fe n] lines have similar shape and brightness, while 
the H2 line is fainter and has a different spatial distribution. 
In particular the emission of [S 11] and [Fe n] lines is maxi- 
mum at knot E with a flux of ~ 4 10~ 14 erg s cirT 2 , while 
the maximum H2 emission is reached in knot P (~ 1.4 10" 14 
erg s _1 cm -2 ). Moreover, the [S n] and [Fen] line intensities 
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Fig. 4. Comparison between the observed 

[Fen]/l7155/[Feii]/l8617 ratios (diamonds) and the theoretical 
curves (dashed lines) computed by Bautista & Pradhan (1996) 
for different values of the electron density, n e , and a fixed 
temperature of 10 4 K. This ratio traces a dense component 
in the jet beam (n e up to 2.5 10 5 cirT 3 in the knots close to 
the source). Note that the zero point of the spatial scale is the 
driving source HH 111 IRS (see Fig|2j. 



HH 111: HIGH DENSITY GAS 





.... T = 5000 K 




y ■ 




T = 10000 K 








T = 15000 K 








nhfjprvAri rnnrjp 


















J 5 * 















4 5 6 7 8 9 

log n e 



Fig. 5. Comparison between observed and predicted 
Caia8540/[Caii]/i7290 ratios. The solid lines indicate 
the range of variation of the observed ratio along the jet. 
The same ratio is computed for three different temperatures 
with a five-level statistical equilibrium code (dotted, dashed 
and dash-dotted curves). From this diagnostic, values of the 
electron density of ~ 6 10 5 - 2 10 6 cm~ 3 are found. 

decrease with distance from the source, while H2 emission is 
stronger in the outer bows (L and P). Such behavior, also noted 
by Nisini et al. (2002), can be explained by the fact that [Fe 11] 
and [S 11] emission is expected to be excited from dissociative 
J-shocks (v s > 30-40 km s _1 ) at the apex of the bows, while the 
H2 emission arises from the bow wings where the transverse 
component of the velocity gives rise to slower C-shocks (v s 
< 30 km s~') that prevent H2 dissociation JDavis et al. 1 996 
IHollenbach 19971 IDavis et al. 199 9 1. Observational confirma- 
tions for such an H2 distribution is seen for example in the HH 
211 flow JGueth et al. 1999} . In this scenario the fact that the 
([Fen] ,[S n] )/H2 observed ratio decreases with distance from 
the source would be consistent with the suggestion that in the 
outer bows there is a prevalence of non-dissociative slow C- 
shocks, since the ambient density becomes higher or compara- 
ble to the jet density JHollenbachl9 89 1. 

In order to apply the diagnostics to line ratios, we first de- 
termined the visual extinction towards all the knots. As ex- 



pected, the extinction is quite low outside the dense core around 
the source, being Ay ~ 2 mag for the innermost visible knot 
(A-B), and decreasing to Ay ~ mag for the farthest knots 
(from II to T-V). These values are much smaller than those de- 
rived in Gredel et al. (1993) and in Nisini et al. (2002), who 
adopted the 1.64/1.25/zm ratio for the Ay determination, and 
smaller but more similar, at least for the inner knots, to the 
values derived from the Balmer decrement (M orse et al. 19931 
pSforiega-Crespo et al. 19 93 1. As explained in Sect. 3.1 and, 
more in detail, in the appendix of Nisini et al. (2005), here we 
assume that the 1.64/1.32/xm ratio yields a more reliable es- 
timate of the reddening than the 1.64/1.25/mi ratio when the 
|Nussbaumer & Storey (198"8")] Einstein coefficients are consid- 
ered. 

After correcting line fluxes for reddening, we derive the 
values of the physical parameters, n e , x e , T e , and n#, sepa- 
rately for the different lines considered, for each knot along the 
jet, with the procedures described in Sect. 3.1. The results are 
shown in Fig. [3] The electron density n e from [S 11] lines (filled 
circles) increases from ~ 7.5 10 2 cm -3 in the first knot, A-B, to 
~ 10 3 cm~ 3 in the brightest optical knots, reaching a maximum 
value of ~ 2 10 3 cirT 3 in F. Then it slowly decreases with dis- 
tance from the source down to n e ~ 0.7 10 2 cm" 3 at knot P, and 
finally presents another local maximum in knot T-V (n e ~ 5 10 2 
cm 4 ). These values agree with the average electron density of 
900 cm 4 found by Hartigan et al. (1994), while they are of a 
factor of two lower with respect to the values found by Morse 
et al. (1993). The electron temperature, T c , derived from the 
optical lines is similar in all the various knots, i.e. cooling re- 
gions along the jet, being on average ~ 1.3 10 4 K, if exception 
is made for the first noisy point at the base of the jet (T e -1.9 
10 4 K). 

The values of the hydrogen ionisation fraction x e , in each 
knot, and averaged over the region of the cooling zones where 
the considered optical lines emit the most, vary from a min- 
imum of 0.05 to 0.17. In our determination, x e decreases 
steadily from knot E to knot Kl, like the electron density, but 
it increases again at knot L to a value of 0.17. In the subse- 
quent faint region of the beam the ionisation fraction appears 
to decrease again, although the accuracy of the derivation is 
much lower here. Finally another local maximum is retrieved 
at the isolated knot T-V very far from the source, where x e is 
0. 14. The lower limit found for the jets brightest knots is sim- 
ilar to the one (x e = 0.052) derived by Hartigan et al. (1994) 
by comparing the observed [N n]/[0 1] ratio, which is sensitive 
to x e , to the predictions of low velocity shock models. In other 
knots, however, we find higher values that are difficult to jus- 
tify, as a shock with effective velocity of 30-40 km s _1 hardly 
produces an ionisation fraction greater than a few percent (see 
Fig.l in Harti gan et al. 1994) . The contradiction can be recon- 
ciled if one considers that the shock fronts in the jet beam ad- 
vance in the wake of the previous shock, and thus move in a 
medium that has been already pre-ionised. We note that a sub- 
stantial level of hydrogen ionisation can persist between two 
consecutive shocks because the recombination time is slow in 
the rarefied jet gas. It can be shown that with the typical elec- 
tron densities and bulk velocities of the jet material, the recom- 
bination time is of the same order of the crossing time of the 
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entire bright section of the jet (see BE99 1. On the other hand, 
this explanation may not apply to the relatively high ionisation 
of knots L and T-V, located at the end of the two sections of 
the jet. In this case a higher ionisation may be produced be- 
cause the shock propagates into a medium of lower pre-shock 
density (see Harti gan et al. 1 994 1. 

From n f and x e we obtain a gross value of the total den- 
sity of 10 3 - 1.9 10 4 ctrT 3 . This is maximum for the brightest 
knots, as expected (the emission in collisionally excited lines is 
proportional to n c -n#). 

The empty circles in the diagrams of Fig. [5] are the values 
of n e and T f inferred from the [Fe n] lines. From these lines we 
find a higher electron density, n e ([Fen]) ~ 1.3 - 4.3 10 3 cirT 3 , 
and a lower temperature, T c ([Fen]) ~ 4800 - 8500 K than in- 
ferred from the optical line ratios. Our result confirms that the 
observed [Fe n] emission typically comes from regions of the 
post-shock zone that are more distant from the front than the 
optical lines, where the gas is cooler and more compressed. 
Such a situation is described in detail in Nisini et al. (2005). 
Note also that our values of n e are lower than those found in 
Nisini et al. (2002) (n, -0.5-2 10 4 cnr 3 ). This is due to the 
fact that here we use the [Fen]1.64yum/1.32/im ratio to deter- 
mine extinction, leading us to infer lower values of Ay- 
Even higher density components can be traced through the 
[Fen]/l7155/[Fen],18617 and Caii/l8540/[Can]/l7290 ratios, as 
illustrated in Sect. 3.1 (Figs.|^and|5|l- From the [Fen] ratio we 
infer values of n e between ~ 10 4 cm -3 and ~ 2.5 10 5 crrT 3 in 
the closest knots E, F, G2, and Gl, that subsequently decreases 
to ~ 10 3 ctrT 3 in the outer knots H, II, J, and L. From the 
[Call] ratio we find n e ~ 6 10 5 - 2 10 6 cm 3 in knots E, F, G2, 
Gl, H, J, while in the other knots no estimate was possible, be- 
cause the signal-to-noise is too low. The values of n e inferred 
from [Fen]A7155/[Fen]/l8617 and CaiLl8540/[Can],17290 ra- 
tios demonstrate that a component of material even denser than 
that traced by [Fen] lines at 1.64, 1.60, and 1.53/mi is present 
in the knots along the jet. If we consider for this higher density 
component the same ionisation fraction found from the optical 
lines, one obtains that n# may be as high as ~ 10 7 cirT 3 in this 
layer of the post-shocked regions. Similar high densities have 
been found in the HH 1 jet JNisini et al. 200 5 ) and at the base 
of a few jets from optical T Tauri stars l |Hartigan et al. 2 004). 

Using the H2 lines flux we derived the conditions of the 
molecular gas in the jet. In the H2 emitting regions the temper- 
ature varies between 2000 and 3000 K, and the column density, 
N//,, between 5 10 16 and 5 10 17 cm~ 2 . From the latter value 
divided by the knot size transverse to the line of sight, we de- 
rive an approximate H2 density of ~ 10 2 cm" 3 . This is only 
a lower limit because it is supposed that the H2 lines are ex- 
cited in a thin layer of molecular material located in the lateral 
wings of the bow shocks (see above). Thus the dimension of 
the emitting region is probably smaller than the knot size. Such 
a picture could be confirmed by high angular resolution im- 
ages. Nevertheless, the H2 lines trace colder and less dense gas 
and for this reason the H2 emission is maximum in the outer 
knots L and P (knot T-V is out of the frame in the NIR range) 
as shown by the spatial profiles of Fig. [2] 
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Fig. 6. Same as Fig. |5] but for the HH 34 jet. The zero point 
of the spatial scale is the source HH 34 IRS (RA(2000): 05 35 
29.8, Dec(2000): -06 26 57). 

The spectacular HH 34 jet is located in the L1641 cloud in 
Orion, and is one of the best studied examples of stellar jets. It 
consists of a parsec-scale flow ( Bally & Devine 1994[ > the first 
30" of which on the blue-shifted side is a well-aligned chain 
of knots emitting in both optical and NIR lines, placed south 
of the central source, HH 34 IRS (Kisloffel & Mundt 1992i . 
HST images of the flow reveal clearly that each knot has the 
morphology of a mini-bow shock (or 'working surface', see, 
e.g., |Ray et al. 199"6||Reipurth et al. 2002| l. Two larger clumpy 
bow-shocks (HH 34N and HH 34S) symmetrically placed at a 
distance of ~ 100" from the source are also visible in the op- 
tical (see, e.g., IMorse et al. 1992L while the counter-jet is not 
detected at optical or infrared wavelengths. Proper motions of 
the knots have been measured from the ground by Eisloffel & 
Mundt (1992) and Pevine et al. (19971 while Reipurth et al. 
(2002) made similar measurements with HST deriving typi- 
cal shock velocities for the working surfaces along the beam 
of only 20 km s . Investigations of the physical properties of 
the jet gas through spectral diagnostic, limited to optical lines, 
have been carried out by Hartigan et al. (1994) and BE99I In 
the infrared band, the jet has been investigated by Stapelfeldt 
et al. (1991), who imaged the jet in [Fen] lines at moderate 
and high angular resolution. The only structure previously re- 
ported north of the star in the H2 2.12 fim line, using a narrow- 
band filter, is an extended arc just south of the optical HH 34N 
JStanke et al.7998l 

We apply our combined diagnostics to the innermost 30" of 
the jet beam where bright optical/NIR lines are seen. We define 
12 knots along the jet, following the optical nomenclature of 
Eisloffel & Mundt (1992). In Fig. [6] we present the emission 
profiles in the optical [S 11] line, the infrared [Fe 11] and H2 
lines, each normalized to its intensity peak. The atomic lines 
are similar both in spatial profile and brightness. The lines in- 
tensity is maximum in knot E (flux ~ 2 10~ 14 erg s~' cirT 2 ) and 
then decreases with distance from the source. The H2 line, on 
the contrary, is fainter and has a different distribution with two 
local maxima in knot C (flux ~ 2.5 10~ 15 erg s _1 cirT 2 ) and 
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Fig. 7. Same as Fig. [3] but for the HH 34 jet. Open circles are 
the values derived from [Fe n] line ratios, while filled circles 
indicate the values determined using the BE technique. Note 
that the zero point of the spatial scale is the driving source HH 
34 IRS (see FiggJ. 



knot E (flux ~ 1.3 10 erg s _1 cm" 2 ). As mentioned for HH 
111, this behavior is due to the fact that [S n] and [Fen] lines 
are excited by the J-shock at the apex of the bows, while H2 
lines are prevalent in slow, non-dissociative shocks, so the rel- 
ative spatial distribution of the lines depend on the prevalence 
of C or J shocks along the jet. We note also faint H2 emission 
between 20" and 31", and between 69" and 80" north of the 
star. The latter region corresponds to the arc detected by Stanke 
et al. (1998). This may indicate the presence of a counterjet that 
is very extincted with respect to the blue-shifted lobe. 



HH 34: HIGH DENSITY GAS 

BCD E FGHIJ K 





3.0 






2.5 




to 






DO 
i< 


2.0 




\ 

in 






m 


1.5 




r-» 






f< 


1.0 




~o 






L_ 


0.5 






0.0 





4 10 s cm"' 



d 1 



2.5_ip'_cm"*_ 
J.6_10'_cm" 5 _ 
10 s cm" 3 



- 0.16 10' cm- 



5 10 15 20 25 30 

distonce from source (arcsec) 

Fig. 8. Same as Fig.|4] but for the HH 34 jet. Values of n e up 
to 4 10 6 cm" 3 are found. Note that the zero point of the spatial 
scale is the driving source HH 34 IRS (see Fig|6j. 
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Fig. 9. Same as Fig. [5] but for the HH 34 jet. The horizontal 
dashed and solid lines indicate the observed range in the knots 
close to the source (A, B, C, D) and in the outer knots (from 
E to L), respectively. Electron densities up to 6 10 6 cm" 3 are 
found in the knots closer to the source. 



The visual extinction (see Sect. 3. 1) is more or less constant 
in the blue lobe with a value Ay = 1.3 mag in the first 30" from 
the source, except for the nearest knot (A), for which we de- 
rive Ay ~ 7.1 mag in agreement with the high extinction value 
found at the source by Reipurth et al. (1986). 

The physical parameters obtained applying the BE tech- 
nique to the optical lines emitted by each cooling zone traced 
by the knots are shown in Fig. Q The inferred values of the 
electron density averaged over each cooling zone, n e , decrease 
from a maximum of 3.4 10 3 cm" 3 in the first knot (A) down 
to 4.3 10 2 cm" 3 in knot D. Then n e increases again up to a lo- 
cal maximum of ~ 10 3 cm" 3 in knot G, and from this point 
on n e starts to decrease again down to ~ 0.5 - 1 10 2 cm" 3 in 
the outer knots K and L. The ionisation fraction, x e , is gener- 
ally quite low in all the knots. It increases from 0.1 1 in the first 
knot, A, up to 0.27 in knot B; then decreases gently all along 
the bright beam down to 0.03 in knots I and J. From there it 
shows isolated higher values up to 0.08 in knots K and L. The 
temperature, T e , derived from the BE technique, turns out to be 
on average ~ 1.4 10 4 K. From x e and n e we obtain an estimate 
of the total hydrogen density, that shows a trend similar to the 
electron density. The maximum n# value is found in knot A 
(n# ~ 3 10 4 cm" 3 ), then increases from 3 10 3 cm" 3 in knot B 
up to 2.7 10 4 cm" 3 in knot G. Subsequently it decreases again 
down to ~ 10 3 cm" 3 in knots K and L. 
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As noticed for the HH 1 1 1 jet in Sect. 3.2, the inferred x e 
and T e values are higher with respect to those obtained for 
the same object in BE99] due to the differences between the 
adopted abundances sets. Nevertheless, the trends are the same. 
The values of n e and x e found for HH 34 agree with the electron 
density derived by Morse et al. (1993), and with the average 
value of the ionisation (x f = 0.026) estimated in Hartigan et al. 
(1994) from the [N n]/[0 1] ratio. 

In Fig.0we also plot the values of n e and T e inferred from 
the [Fe n] lines. As in the case of HH 111, n e derived from the 
[Fe n] lines is higher than obtained from the [S n] lines {n e ~ 1 

- 4.5 10 3 cm" 3 ), while the values for T e are lower (T e ~ 3500 

- 7700 K). Such behaviour again corresponds to the expected 
structure of the shock cooling zones. In this jet too we find an 
even denser component of gas in the knots, as shown in Figs. [8] 
and[9] From the [Fen]/l7155//18617 ratio we infer values from 
n e 
L. 



~ 4 10 6 ctrT 3 close to the source to n e ~ 10 2 cirT 3 in knot 
Similar results are obtained through Can diagnostics (Fig. 

- 5 - 6.3 10 6 cm -3 for the inner 
10 6 cirT 3 in the outer knots (from 
E to L). Further estimates of relevant physical quantities are 
discussed in Sect. 5. 



In this case we have n e 
knots (from A to D) and n e < 



4.3. HH 83 

The HH 83 jet is also located in the Orion Nebula, in the L1641 
molecular cloud (Reipurth 1989). It is powered by the IRAS 
source 05311-0631, which is detected at infrared wavelengths 
only from the H and K band ( |Moneti and Reipurth 19 95 1. The 
jet becomes visible in the optical and near-infrared at a distance 
of ~ 10" from the source, when it emerges from the cloud sur- 
face, detected as a reflection nebula (Re 17, Rol ph et al. 19 90 1. 
The jet is composed of 10 blue-shifted knots (following the 
nomenclature of Reipurth 1989 1 located between 10" and 50" 
from the source, and terminates with a bow-shock at ~120". 
On the opposite side of the source only two faint condensa- 
tions are found (Reipurth 1989 1. Because of the faintness of 
the H2 emission in the acquisition image, the slit was not well 
aligned for the infrared spectra. This caused the IR spectra to 
be degraded, and only a few, very faint [Fen] and H2 lines 
were observed. Therefore it was not possible to extract the IR 
spectra for combination with the optical ones. We concentrate 
on the optical emission from the blue lobe of the jet, applying 
only the optical diagnostics. From the optical [Sn], [Oi], and 
[Nn] lines, using the BE technique we derive the values of n e , 
x e , T e , and n#. These are shown in Fig. EH an d are plotted ver- 
sus the distance from the driving source. The electron density 
reaches a maximum value of ~ 0.7 10 3 ctrT 3 in knot E; then 
it decreases with distance from the star down to ~ 70 cirT 3 in 
knots I-J. The ionisation fraction and the temperature could not 
be calculated in the first three knots (A-B-C, D, E) because the 
diagnostics did not find a solution for the measured line fluxes. 
This is probably due to the presence of the reflection nebula, 
the continuum from which could not be properly subtracted, 
adding spurious flux to the jet emission lines. Concerning the 
other knots, the ionisation fraction has a maximum value of ~ 
0.55 in knot F, and then it decreases down to ~ 0.22 in knot 
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Fig. 10. Same as Fig. but for the HH 83 jet. No infrared 
diagnostics are available in this case (see text). The zero point 
of the spatial scale is the source HH 83 IRS (RA(2000): 05 33 
32.5, Dec(2000): -06 29 44). 



H. There is another local maximum in knot I-J, for which we 
found x e ~ 0.46. The temperature is relatively high, T e ~ 2.1 
10 4 K in knot F and lower in the following knots, down to 1.4 
10 4 K in knot H. Then it is higher again in knot I-J (T f ~ 1.8 
10 4 K). Finally, the total density, n#, assumes its highest value 
of ~ 10 3 cirT 3 in knots F and G. We note that a relatively higher 
degree of ionisation is found for this and the following jets with 
respect to HH 34 and HH 111. For an interpretation of this re- 
sult, as well as for the determination of the mass flux in this 
flow, see the discussion in Sections 5.2 and 5.3. 
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4.4. HH 73 
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Fig. 11. Same as Fig. [10] but for the HH 73 jet. Due to low 
signal-to-noise, the results are subject to large errors. Since the 
source of this jet is not known, the zero point of the spatial scale 
was set at the beginning of the first detected knot (knot A). 



The HH 73 jet, discovered by |Reipurth & Graham ( 1988) 
is located in the Vela Molecular Ridge (D = 450 pc), near the 
small nebula Re 6 (Reipurth 198 1| . It is composed of a chain 
of very faint optical knots, that we name from A to G, and ex- 
tends for a total length of ~ 70". Since the source of this jet 
is not known, the zero point of the spatial scale was set at the 
beginning of knot A. We do not detect any emission in the near 
infrared, so we applied our diagnostics only to the optical emis- 
sion of the brightest knots (from B to F). The derived values of 



n e , x e , T e , and n# are shown in Fig. Because of the low 
signal to noise of the line fluxes, our diagnostics is affected by 
large errors. Nevertheless, we can define a range of variation 
for the physical parameters. The values of n e are between ~ 0.2 
10 3 cirT 3 and 1.4 10 3 cirT 3 ,x e varies between 0. 1 and 0.45, the 



1.4 10 4 - 2.2 10 4 K, and the 
1-4 10 3 cirT 3 . The relatively 
high ionisation in this jet may be justified following the same 
arguments identified for HH 83 (see Sect. 5.3). 



temperature T e varies between 
total density is in the range of - 



4.5. The HH 24 complex 

The HH 24 complex lies in the NGC 2068 nebula that 
includes both reflection and emission nebulae, and at 
least three different outflows JSolf 19871 IMundtetal. 19911 
Eisl offel & Mundt 19971 IBTMI. IR imaging by |Lane(1989)| 
and |Zealey et al. (1 989) revealed the presence of three opti- 
cally invisible sources (SSV 63, SSV 63 W, and SSV 63 NE) 
which could be the driving sources of the observed outflows. 
One of the best defined structures is the HH24 C blue-shifted 
jet, aligned with the source SSV 63. Eisloffel & Mundt (1997) 
showed that the HH 24 C blue lobe extends for ~1 pc, includ- 
ing the ring-shaped object called HH 24 H, and terminating 
with a giant bow formed by a number of HH objects (HH 20, 
HH 21, HH 37, HH 70). On the opposite side there is a chain of 
red-shifted HH knots (knots of the group E, knot A, and knot 
M following the nomenclature of Eisloffel & Mundt 1997) that 
seems to form the counterjet, although the structure of this re- 
gion is still unclear. The red lobe, in fact, deviates from the HH 
24 C flow axis by ~ 6° and there is evidence that the bright knot 
A may belong to another jet in the complex (HH 24-MMS jet 
according to Eisloffel & Mundt 19971 . In our spectra we de- 
tected the innermost ~80" of the blue lobe (from knot CI to 
knot CI 1) and -30" of the red lobe (knots El, E2, and A). In 
the optical many lines are visible throughout all the jet length, 
while in the infrared we observed only a few faint [Fe n] and 
H2 lines (from the red-shifted knots El and A). As a conse- 
quence, our diagnostics is limited to the derivation of the phys- 
ical parameters in the region of optical emission through the 
BE technique (see Fig. I12i. In the blue lobe the electron den- 
sity n e decreases with distance from the source from ~ 8 10 2 
cm~ 3 down to ~ 2 10 2 cm" 3 in knots C5-C6. Then there is a 
local maximum in knot C7 (n c ~ 4 10 2 cm -3 ). As for the red 
lobe, n e is equal to ~ 7 10 2 cirT 3 in knot El, then it decreases 
down to ~ 3 10 2 cirT 3 in knot E2, and finally it reaches a value 
of ~ 8.5 10 2 cirT 3 in knot A. The ionisation fraction x e varies 
along the jet between 0.2-0.4, with no definite trend. It is lower 
in knot El, where we find x e ~ 0.05. The temperature decreases 
with distance from the source in the blue lobe going from -2.1 
10 4 K down to ~ 1.7 10 4 K. In the red lobe, instead, we found 
a temperature of ~ 2 10 4 K in knot El, ~ 1.4 10 4 K in knot 
E2, and a maximum of ~ 2.4 10 4 K in knot A. Finally, the total 
density is decreasing in the blue lobe from ~ 3 10 3 cirT 3 in knot 
CI to 6 10 2 cm -3 in knot C5, then it increases again in knots C6 
and C7 (ny ~ 0.7-1.2 10 3 cm 3 ). Towards the red lobe we find 
Tin ~ 1.4 10 4 cirT 3 in knot El, n# ~ 10 3 cirT 3 in knot E2, and 
Tin ~ 4 10 3 cirT 3 in knot A. We find higher ionisation fractions 
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Fig. 12. Same as Fig. EH but for the HH 24 C/E jet. The zero 
point of the spatial scale is the source SSV 63, which is sup- 
posed to be the driving source of this jet (RA(2000): 05 37 
08.23, Dec(2000): -00 09 24.62). 

and temperatures than those derived by BE99 As mentioned 
in the case of the HH 34 jet, this is due to the adoption of a 
different set of abundances (see also Fig. 1). For comments on 
the comparison with HH 111 and HH 34 see the discussion in 
Sect. 5.3. 

It can be noticed that, while in the blue lobe the derived 
parameters have a definite trend along the jet, in the red lobe 
the values are scattered. We find, e.g., that the ionisation frac- 
tion and the total density in knot El are, respectively, much 
lower and higher than those found in the C knots. This effect 
can be due to the faintness of knots El and E2 (see the inten- 
sity spatial profiles in Fig. I12> and, as a consequence, to the 



low signal-to-noise ratio, which affects, overall, the measure- 
ment of the [Nn]/16583 line flux and, thus, the determination 
of the ionisation fraction. It could also be intrinsic to the ob- 
ject, the structure of which is not entirely clear in this region 
(Sol fT987l IMundt et al. 19911 lEisloffel & M undt 1997) (knot 
A may not belong to the jet and the emission from the HH 24 
J jet can be superimposed on the knots of group E). See further 
discussion about the HH 24 C/E jet in Sections 5.2 and 5.3. 

Another jet in the complex is HH 24 J. It consists of 
a chain of knots and condensations, which are well aligned 
with the source SSV 63 W. This jet also seems to be 
parsec-scale, terminating with the bow HH 19 to the north- 
west and the bow-shaped object HH 27 to the south-east 
(Eisloffel & Mundt 19971 . In the optical we detect emission in 
the first ~30" of the north-west lobe and in the first ~40" of the 
south-east lobe (HH 24 J and HH 24 K following the nomen- 
clature of Eisloffel & Mundt 19 97 1, while no emission lines are 
detected in the NIR. Unfortunately, in the spectra of the knots, 
the [S ii]/1/16716, 6731 doublet is not resolved sufficiently as to 
separately measure the line fluxes, except for two condensa- 
tions located at ~30" on the north-west from the source. Since 
higher resolution spectra of the HH 24 J jet are not available, 
we apply our diagnostics only to the two condensations. For 
these knots we derive an electron density of ~0.5 10 2 cirT 3 , a 
ionisation fraction of ~0.02 and -0.05, a temperature of ~1.5 
10 4 K and -0.9 10 4 K, and a total density of -2.5 10 3 cm 3 and 
0.9 10 3 cm" 3 . 



5. Discussion 

Here we illustrate and discuss further information about the 
jet physics, derived from the basic gas parameters obtained in 
Sect. 4 and summarized in Tab. [5] We also compare the values 
of these parameters in the various objects, to search for com- 
mon trends or differences in the sample. 

5.1. Abundances of refractory species 

The wavelength range covered by our spectra includes tran- 
sitions from refractory species like carbon, calcium and iron. 
This allows us to check for the gas phase abundances of these 
species with respect to solar values. In the interstellar medium 
refractory species are often locked into dust grains, thus a 
strong depletion is expected. On the other hand, sputtering and 
photoevaporation processes due to the passage of shock fronts 
can destroy all or part of dust grains, releasing the refractory 
atoms into the gas cloud dJones 2 000 Drai ne 2003> . It follows 
that an estimate of the depletion of the gas phase abundance of 
these atoms can give important constraints on the dust struc- 
ture, and on the efficiency of shocks in destroying grains. The 
latter is related to parameters such as the shock velocity and the 
pre-shock excitation conditions. In order to infer the elemen- 
tal gas-phase abundances we compare observed and expected 
ratios between emission lines of refractory and non-refractory 
species. We select lines which are excited in the same region of 
the post-shocked gas, so that one can assume the same filling 
factor. 
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Table 2. Physical parameters of the jets averaged over the brightest knots. 
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(10 3 K) 


HH 111 


2.0-0 


1.0 


0.10 


13.0 


11.3 


2.6 


7.3 


HH 34 


7.1 - 1.3 


0.8 


0.04 


13.8 


16.2 


1.8 


5.8 


HH 83 




0.5 


0.38 


17.5 


0.9 






HH 73 




0.6 


0.31 


17.5 


1.7 






HH 24 C/E 




0.4 


0.32 


19.3 
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Table 3. Gas-phase abundance of refractory species with respect to the Solar Abundances determined by Asplund et al. 2005 



Species 




[X] sas /[X] 


saUtr 


Abundances: (X/H) 




HH 111 


HH 34 


Orion Cloud 


Orion Abundances" 


Solar Abundances* 


Ca 


0.5 - 1 


0.3- 1 


0.01 


2 10~ s 


2.04 10~ 6 


C 


~ 1 


~ 1 


1.07 


2.63 10~ 4 


2.45 10" 4 


Fe 




0.13 


0.05 


1.29 10~ 6 


2.82 10~ 5 







"from Esteban et al. (2004), except for the Ca abundance which was determined by Baldwin et al. (1991) 
''from Asplund et al. (2005) 



For the gas phase abundances of calcium and carbon 
we use the ratios [Can]A(7290+7324)/[Sn]/l(6716+6731) and 
[Ci]/l(9824+9850)/[Sii]/l(6716+6731). The predicted values 
of these ratios are determined using the values of n e , x c , 
and T e inferred through the BE technique. Calcium is as- 
sumed to be completely ionised while the carbon ion- 
isation fraction is computed considering collisional and 
charge exchange ionisation and direct and dielectronic re- 
combination processes (rates from IStancil et al. 1 998 and 
|Landini & Monsignori Fossi 1990) . Finally, we assumed solar 
abundances from Asplund et al. (2005) for Ca and C. In Fig. 1131 
the comparison between predicted and observed ratios along 
the HH 1 1 1 and HH 34 jets is shown. 

The observed [Caii]/l(7290+7324)/[SiiR6716+6731) ra- 
tio is smaller than the predicted one along both HH 34 and 
HH 111. We interpret this discrepancy as due to depletion of 
Calcium atoms with respect to the solar value (see Tab. [3}. 
Along the HH 1 1 1 jet a depletion of ~ 50% is found in the knots 
closest to the source. Then the depletion slowly decreases with 
distance, down to zero in knot O. Another local maximum is 
then found in knot P. Similarly, along HH 34 we find a maxi- 
mum depletion (~70%) in the first knot, which then decreases 
steadily until a solar gas-phase abundance is reached in the out- 
ermost three knots (J, K, L). 

The trend, observed along both HH 111 and HH 34, can 
be explained by the fact that the leading bow-shock has a 
much larger effective shock velocity (of 100-400 km s~') than 
the internal working surfaces in the beam (30-40 km s~') 
(Raga& Kofman 1992 1. On the other hand, the shock mod- 



els predict substantial grain destruction only for shock speeds 
> 100 km s _1 dDraine 20 03). We would therefore only expect 
the dust to be gradually destroyed by the working surface as it 
propagate outwards. At leading bow-shocks however, e.g. HH 
34 S, almost total destruction should be achieved. Note that 
along the jet we measure less depletion than in the surround- 
ing ambient medium (see, e.g., Bald win et al. 19 91 1, which is 
another indication that dust is partially destroyed by jet shocks. 

For carbon, instead, good agreement between the observed 
and predicted [C i]/l(9824+9850)/[Sn]/l(6716+6731) ratio is 
found, at least in the brightest knots where the signal-to- 
noise is high. This indicates that this species is not depleted 
in the Orion Cloud, a result that has also been found in 
previous works (e.g., IRubin et al. 19911 Baldwin et al. 199TI 
Pei mbert et al. 1993lEsteban et al. 20041 . 

In a similar way we determine the quantity of Iron 
atoms in gaseous form. For Fe the comparison is compli- 
cated by the fact that the observed [Fen] lines come from 
regions of the post-shocked gas different from those giving 
rise to the optical lines. This is why contrasting results are 
found in the literature, in works that adopt different meth- 
ods (see e.g. IBeck-Winchatz et al. 19961 IBohmetal. 20011 
Nisi ni et al. 20 02 1. Since we obtained different values of n e 
and T e from [Fe n] and [S n] lines, we cannot compare ra- 
tios between these lines, because they have probably different 
filling factors (see N isini et al. 20051 . We use, instead, the ra- 
tio [Feii]1.25/im/[Pn]1.18/im. Phosphorus, in fact, is a non- 
refractory species, and the excitation conditions of the 1.18/mi 
line are similar to the ones of the [Fen]1.25^/m line. The 
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Fig. 14. Comparison between the observed 

[Fen]1.25//m/[Pn]1.18/im ratios (crosses) and the ex- 
pected value obtained assuming solar abundances for Fe and 
P JOliva et al. 2001 1 . The position of the observed points 
indicates a depletion of ~ 85 % of the gas phase abundance of 
Iron. 
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Fig. 13. Comparison between observed (dia- 

monds) and predicted (dashed and dash-dotted 
lines) [Can]/l(7290+7324)/[Sii]^(6716+6731) and 
[Ci]^(9824+9850)/[S n] ^(67 16+6731) ratios along the 
HH 111 and HH 34 jets. The theoretical ratios are computed 
from the derived physical parameters of the gas and assuming 
solar abundances for Ca and C. The [Ca n]/[S n] ratios indicate 
a depletion of the gas phase abundance of Ca of ~ 50% in 
the innermost knots. The depletion decreases with distance 
from the source for both HH 111 and HH 34. The good 
agreement between the observed and theoretical [Ci]/[Sii] 
ratios indicates that carbon is not depleted with respect to 
solar values. Note the vertical broken lines in the HH 111 plot 
indicate a gap in the distance axis. 



predicted [Fen]1.25yum/[Pn]1.18jum is estimated to be about 
[(Fe/H)/(P/H)]/2 dOlivaetal. 200 » if Iron is not depleted. The 
comparison between the observed ratio in the HH 34 jet and the 
one calculated assuming solar abundances is shown in Fig. 1141 
From the figure we deduce that in HH 34 Fe may be depleted 
by 87% with respect to solar ((Fe/H) faj /(Fe/H) JO ; ar ~ 0.13 and 



thus (Fe/H)g fls ~ 3.67 10~ 6 ). This depletion is lower than the 
one found by Esteban et al. (2004) in the ISM of Orion, so 
some atoms of Iron have probably been unlocked by the jet 
shocks. The fact that we found a different amount of deple- 
tion for Calcium and Iron indicates that the different species 
follow selective patterns for their erosion from dust grains as 
expected from theory (Jones 2000 1. Similar results have been 
also obtained on the HH 1 jet in Nisini et al. (2005) indicating 
a common behavior of the erosion patterns in different jets. We 
could not check the iron depletion along the HH 1 1 1 jet, since 
the [P n] 1 . 1 8yum line is not visible in our spectra. 

5.2. Determination of the mass flux rate and other jet 
parameters 

The mass flux rate (Mj e t) is a fundamental quantity gov- 
erning the jet dynamics and which enters all comparisons 
between observations and theoretical models. For example, 
in the magneto-hydro dynamic models proposed to explain 
jet formation and acceleration (see [Konigl & Pudritz 2 000 
IShu et al. 2000> . the ratio between the rate of mass ejected 
into the jet (Mj et ) and the rate of mass accreted from the 
disk onto the star (M acc ) is fixed (Mj et /M acc ~ 0.01 - 0.1). 
Moreover, the knowledge of Mj et allows us to estimate other 
important dynamical quantities such as the linear (Pjet) and 
angular (Lj et ) momentum fluxes carried by the jet. Knowing 
Pjet one can check if the jet is powerful enough to accelerate 
surrounding molecular outflows, thus helping to clear the cir- 
cumstellar environment and inject turbulence into the cloud. 
Lj et is related instead to the jet's capability of removing ex- 
cess angular momentum from the disk/star system, thus allow- 
ing the accretion of matter from the disk onto the central star 
( |Konigl & Pudritz 2000|IWoitas et al. 20051 . 

The determination of the physical conditions along the jet 
allows us to derive Mj et in each knot. We use two different pro- 
cedures, as we also did in Nisini et al. (2005) for the HH 1 
jet. In the first method (method A), we estimate the mass flux 
as Mj et = HtnHiiHitrj 1 vj, where /U=1.24 is the mean atomic 
weight, rrin the proton mass, n# the hydrogen density and rj 
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Table 4. Mass flux and linear momentum flux along the jets 
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IMundt et al. T99T1 (HH 24 C). 

V jet is taken to be one half of the FWHM of the [Sn] intensity profile measured by |Reipurth et al. 2000"| (HH 1 1 l), |Reipurth et al. 2002| (HH 
34), IMundt et al. 1991K HH 83, HH 24 C). 

c n H is the total density estimated via the BE technique (see Sect. 3.2). 

d ff is the volume filling factor estimated from the ratio between observed and theoretical line luminosities; this is equal to the ratio of Mj et 
values derived from methods B and A (see text). 

f Mj5. t estimated through method A (see text) and averaged over the brightest knots; a //=1 is assumed. 

^Mj ct measured from [S n] and [Oi] line luminosities taking into account the beam filling (method B) and averaged over the brightest knots. 



:M iE 



s Pj ct is calculated from Mj ct (B) and vj ct (.Pj ct 
''^outflow is the flux of linear momentum tr 
Reip urth & Olberg 1991||Cernicharo & Reipurth 199g| (HH 111). 



''^outflow is the flux of linear momentum transported by the molecular CO flows and measured by IChernin &Masson 19 95 (HH 34) and 



and v y, respectively, the jet radius and velocity. We use the total 
density inferred from our diagnostics and the values of rj and 
v j taken from the literature, as our spectral resolution is too low 
for such an estimate (rj is taken to be one half of the FWHM 
of the [S n] intensity profile in HST images). This method is 
independent of the reddening estimate, but it assumes that the 
knot is uniformly filled at the density derived from the diag- 
nostic, giving an upper limit to Mj et . Such an effect is partially 
compensated for by the presence of regions at even higher den- 
sities in the beams than those traced by the [S n] lines, which 
are not taken into account in the calculations. Alternatively 
(method B), Mj et can be derived from the observed luminos- 
ity, L(line), of selected optically thin lines such as [Sn], [Oi], 
and [Fe n] , that is proportional to the mass of the emitting gas: 

M jet = pL m H (n H V) v,ll u with n H V = L(line) (li v A; /iff) , 
where V is the volume effectively filled by the emitting gas, 
v, and I, are the tangential velocity and length of the knot, A, 
and fi the radiative rate and the upper level population relative 
to the considered transition, y and f are the ionisation frac- 
tion and the relative abundance of the considered species. This 
method is affected by uncertainties in absolute calibrations, ex- 
tinction, and distance, but does implicitly take into account the 
volume filling factor //, which, in practice, is simply the ra- 
tio of Mjet values derived from methods B and A. In fact, the 
method only measures gas which is sufficiently heated to radi- 
ate the observed lines. 

The mass fluxes inferred for the HH 111, HH 34, and 
HH 83 jets are shown in Fig. [B] For HH 111, method A 
gives a mass flux approximately constant with distance from 
the source and equal to about ~ 2 1(T 7 M yr _1 (rj are from 
|Reipurth et al. 2000| and v, from |Hartigan et al. 200T) l. Method 
B gives smaller values of Mj et , that, however, decrease with dis- 
tance from the source. A steady jet flow is actually expected to 



show a more or less constant Mj et along the beam (as derived 
from method A) unless significant mass is loss sideways by, 
e.g., a turbulent boundary layer or big bow shocks. This seems 
to be the case for HH 1 1 1 . A possible explanation for the ob- 
served decrease is that the jet opening angle increases with dis- 
tance from the source to the point that the jet diameter becomes 
larger than the slit width from 30" outwards. Thus method B 
determinations are affected by partial loss of line flux. To ob- 
tain a value of the mass flux which takes into account the beam 
filling, we use the filling factor from the knots in the first 30" 
(ff~ 0-2) to correct the values obtained from method A. We 
obtain an average Mj et value of ~ 5 10~ 8 M yr~'. Finally, we 
estimate the mass flux transported by the molecular component 
through the inferred H2 column density and the velocity of the 
molecular gas measured by Davis et al. (2001). From method 
A we obtain Mje t (H 2 ) ~ 10~ 9 -10~ 10 M yr _1 , that is down two 
orders of magnitude with respect to the atomic component. As 
explained in the previous Section, H2 emission is thought to 
arise in the C-shocks that form in the lateral wings of the bow 
shocks. Although high angular resolution images in these lines 
are not yet available to confirm this picture, the low value we 
find for the mass loss rate is consistent with H2 tracing only a 
thin shocked layer of the jet. 

For HH 34 we obtain again different values from methods 
A and B, (rj and v, are from Reipur th et al. 200 2 ) but in this 
case both are constant with distance, as the jet beam is always 
smaller than the slit width. Note that the values of Mj et derived 
from [Fen] luminosity are a lower limit since in this calcula- 
tion we assumed that all Fe is in gaseous form. Taking into ac- 
count the beam filling, and considering only the brightest knots 
(E-J) for which we have good signal-to-noise, we infer an av- 
erage value of the mass flux of ~ 3.9 10~ 8 M yr~'. It can be 
noticed that both for HH 1 1 1 and HH 34 the values of M Kt de- 
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Fig. 15. Mass fluxes determined for the HH 111, HH 34, and 
HH 83 jets. Triangles refer to values of Mj et inferred from our 
estimates of the total density, n#, and the values of the jet ra- 
dius (rj) and velocity (vj) taken from the literature (method 
A). Circles (S + , N + , O ) and crosses ([Fen] ) refer to the esti- 
mates of Mjet derived from line luminosities (method B). This 
method takes into account the knot beam filling. Squares re- 
fer to the values of the mass flux transported by the molecular 
component. These values are calculated using our estimates of 
H2 column density (method A). 



rived from method A are similar to the ones found by Hartigan 
et al. (1994) using the same technique. The values obtained 
taking into account the volume filling factor (method B) are in 
both cases lower by one order of magnitude. The latter values 
are a better estimate of the mass flux transported by the atomic 
optical component. Nevertheless, they are lower limits to the 
total mass flux since we are not considering even denser com- 
ponents nor the material that does not radiate (as shown, e.g., 
by the Mj et ([Feii] ) estimates). 



For HH 83 we infer the mass flux from the derived total 
density and the values of the jet radius and velocity from Mundt 
et al. (1991) and Reipurth (1989). We could not estimate the 
mass flux from line luminosities because we do not have a mea- 
sure of the visual extinction, due to the non-detection of NIR 
lines. The mass flux from method A is Mj et ~ 6.9 1CT 8 M Q yr _1 . 
Lacking an estimate of // here, this is only an upper limit. 

For HH 73, and HH 24 J, no literature determinations of the 
jet radius, velocity, and the angle of inclination with respect to 
the plane of the sky are available. Thus we could not estimate 
the mass flux for these objects. 

Finally, for HH 24 C/E we obtain raw estimates of the av- 
erage mass flux from the average jet radius and radial velocity 
given by Mundt et al. (1991). Combining our average value of 
the total density and assuming an inclination angle with respect 
to the plane of the sky of 25° (BE99l we find a mass flux of ~ 
9.9 lO^Moyr 1 . 

For our sample of HH jets we determine the linear momen- 
tum as Pjet = Mjet Vjet, combining the estimated Mjet with the 
average values of vj taken from the literature. We find Pj e t ~ 
1.3 10" 5 Moyr-'kms- 1 for HH 111, P jel ~ 0.8 1(T 5 Moyr" 1 
km s" 1 for HH 34, Pj e , ~ 1.5 1(T 5 M^r 1 km s 1 for HH 83, 
and Pjet ~ 4.2 1(T 5 Moyr 1 km s" 1 for HH 24 C. The results 
are summarized in Tab. |4] 

For HH 111 and HH 34 we compare these values with those 
measured for the molecular outflows seen in CO lines. These 
were estimated by Chernin & Masson (1995) for HH 34 and 
Reipurth & Olberg (1991) and Cernicharo & Reipurth (1996) 
for HH 111. The comparison (see Tab. |4} shows that the flux 
of linear momentum carried by the HH 111 and HH 34 jets is 
higher or comparable to that of the molecular outflows. Thus 
in principle the jets appear to be capable of accelerating these 
outflows. 

In order to measure the flux of angular momentum 
carried away by the jet, we need measurements of the 
jet toroidal velocity, which requires sub-arcsecond resolu- 
tion spectra (see, e.g,. IBacciotti et al. 20021 |Coffey et al~2 004 
IWoitas et al. 2005V To this aim one should acquire high angu- 
lar and spectral resolution spectra with ground-based instru- 
ments equipped with adaptive optics. 



5.3. Comparison between the examined objects 

Our analysis allows us to compare physical parameters of the 
jets in our sample. In particular, it is interesting to investigate 
whether the various excitation conditions in our sample arise 
from different characteristics intrinsic to the jets, or are set by 
the environment through which the jets propagate. From the 
analysis of the reduced spectra we can divide the examined tar- 
gets into two classes: (I) jets visible both in the optical and in 
the near-infrared (i.e. HH 1 1 1 and HH 34) and (II) jets that 
show few, very faint or no lines in the near-infrared (i.e. HH 
83, HH 73 and HH 24 C/E). The parameters collected in Tab.|2] 
show that the jets which are not visible or faint in the infrared 
are less dense than those that show NIR emission. This is ex- 
pected, since the Fe + levels from which infrared lines originate, 
have critical densities >10 4 cm 4 . Thus, electron densities >10 3 
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cm" 3 are needed to populate them. Jets in group (II) also show 
higher ionisation fractions and temperatures. The fact that the 
ionisation fractions are higher in the less dense jets is consistent 
with the shocks propagating in a medium of low preshock den- 
sity. This produces a higher excitation for a given shock veloc- 
ity plartigan et al. 1994) . The low electron density also slows 
down the recombination process, keeping the level of ionisa- 
tion higher. An external source of UV radiation could also con- 
tribute to the excitation conditions. Such a source, however, re- 
mains unidentified. 

The above arguments would indicate that the ionizing agent 
is identical and has the same efficiency, in both groups of jets. 
On the other hand, the lower temperatures inferred for Group 
(I) jets indicate that the cooling process is more efficient in the 
denser objects, as expected for collisional line excitation. 

We also note that the mass flux does not vary significantly 
for the two samples. If we expect a steap decrease of the mass 
flux with age, it seems that the sources of the two different 
groups have similar ages. 

HH 111 and HH 34 are two examples of giant HH flows 
and have similar characteristics such as: (i) a well defined 
blue lobe consisting of a long chain of bright knots and an 
almost completely obscured red lobe; (ii) the presence of dis- 
tant bow-shocks, such as the bow T-V along the HH 111 blue 
lobe at a distance of 150" from the source and the HH 34 
N a nd HH 34 S bows at a distance of 100" from HH 34 
IRS (Reipurth & Bally 2001 1); (iii) similar velocities (see Tab. 
|4j; (iv) similar [Sn] and [Fen] profiles decreasing with dis- 
tance from the source; (v) similar physical gas conditions (see 
Tab. (vi) both jets have a denser component traced by 
[Can] and [Fen] lines; (vii) along both jets Carbon is of so- 
lar abundance, while Calcium is depleted. Moreover the de- 
pletion has a similar trend, decreasing with the distance from 
the source; (viii) similar flux of mass and linear momentum 
(see Tab. 0}. On the other hand, they show important differ- 
ences. The HH 34 jet is less extincted and both the source 
and the inner knots are visible in the optical, while the HH 
111 source is deeply embedded in the parental cloud and the 
jet becomes visible only at a distance of ~ 15" from the 
source. As a consequence the HH 34 jet is associated with 
a small molecular outflow, while along the HH 111 axis a 
powerful CO flow has been detected dChernin & Masson 19951 
|Cernicharo & Reipurth 199~6) . Moreover, whereas HH 1 1 1 is a 
strong H2 emitter, with a maximum H2 2.12 /im line flux of ~ 
1.4 10 erg s _1 cirT 2 from knots L and P, the HH 34 jet shows 
only a few weak H2 lines, with a maximum flux of 2.5 10 -15 
erg s _1 cirT 2 along the 2.12 /jm line from knot C. This differ- 
ence is probably due to the fact that the environment surround- 
ing HH 34 is less dense than the one surrounding HH 111. In 
fact, H2 lines are likely to be excited in the bow wings, where 
the ambient density encountered by the jet is higher than, or 
comparable to, the jet density and low-velocity shocks (v ~ 20 
km s" 1 ) are prevalent, preventing H2 dissociation. 

6. Summary and conclusions 

The application of our combined optical/NIR diagnostics to a 
sample of protostellar jets allows us to define their physical 



structure, to infer important parameters related to their dynam- 
ics such as their mass and linear momentum fluxes, and to in- 
vestigate the acceleration of associated molecular outflows, as 
well as the mechanism of dust reprocessing. For the HH 111 
and HH 34 jets, which are detected both in the optical and in 
the near infrared, we obtain the following results: 

- In HH 34 the visual extinction, Ay, derived from [Fen] 
lines varies between Ay~7.1 mag, at a distance of 3" from 
the source, to 1.3 mag, from 5" to 30" from the source. 
Along HH 111 the visual extinction is Ay~2 mag between 
15"-30" from the source and negligible in the outer knots. 

- In both jets the electron density from the [S n] lines varies 
between 0.5-3.4 10 3 cirT 3 , while in the more compressed 
region traced by [Fe n] lines we find higher electron den- 
sity (n e ~ 1-5 10 3 cm 4 ). An even denser layer with electron 
densities up to 10 6 cm -3 is traced by [Can] lines. Such a 
density stratification is expected in the post-shocked gas of 
each unresolved knot. 

- A temperature stratification is also found: the temperature 
derived from optical S + , N + and O lines is on average 1.3 - 
1.4 10 4 K, while in the region of Iron emission the temper- 
ature is always lower than 10 4 K. 

- In HH 111, which is a strong H2 emitter, we also derive 
the physical conditions of the molecular component. These 
lines probably trace the lateral wings of the bow shocks 
associated with the knot working surfaces, where H2 post- 
shock temperatures of 2000 - 3000 K are found. 

- The ionisation fraction derived through the BE technique 
demonstrates that the gas in the jet is only partially ionised 
(x c ~ 0.03~0.3). The derived total density (n e /x c ) is about 
0.1-3 10 4 cm" 3 in the region of optical emission. If we as- 
sume the same ionisation fraction in the more compressed 
regions where [Fe 11] and [Ca 11] lines are excited, total den- 
sities up to 10 7 cm 3 are obtained in the dense layers traced 
by these lines. 

- Estimates of the gas-phase abundance of carbon, calcium 
and iron are obtained. We find that carbon is not depleted 
with respect to solar values, while calcium and iron are 
strongly depleted (the depletion is ~ 87% for iron atoms 
and varies between 70% and 0% for calcium). On the other 
hand, the estimated depletion is lower than the one found 
in the Orion Cloud by Baldwin et al. (1991) and Esteban 
et al. (2004) both for Calcium and Iron. This result demon- 
strates that weak shocks only partially destroy dust in the jet 
and that some grains still survive. In fact a velocity of 100- 
400 km s _1 , like that of a primary bow-shock, is required to 
completely destroy the dust ( Draine 20Q3jl. 

- Taking into account the filling factor of the gas in the knots, 
we derive mass flux rates of ~ 4-5 1 0~ 8 M yr" 1 . In HH 1 1 1 
the flux of mass transported by the molecular H2 compo- 
nent is two orders of magnitude lower than the atomic one. 

- The estimated flux of linear momentum is on average ~ 
10~ 5 M yr _1 kms _1 . This is higher than, or comparable 
to, the flux of momentum measured for the associated CO 
molecular outflows, and indicates that jets can drive these 
molecular outflows. 
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For the HH 83, HH 73 and HH 24 C/E jets we could only 
apply optical diagnostics, obtaining the following results: 

- The electron densities inferred from the [S n] doublet are 
very low (n e < 10 3 cm 3 ). This explains why these HH jets 
are not detected in the NIR. In fact, the transition of [Fe n] 
lines comes from levels which have a critical density higher 
than 10 4 cm 4 . 

- The ionisation fraction, on the other hand, is always higher 
in these jets (x e up to 0.6), consistent with the fact that for a 
given shock velocity a higher ionisation is produced when 
the shock propagates into a low density medium, and with 
the fact that recombination is slowed down at low electron 
densities. The latter can also be partially responsible for the 
higher temperatures in these objects (T c ~ 1.4-3 10 4 K), 
because the collisional cooling process is less efficient at 
low density. We suggest that the efficiency of the excitation 
mechanism is the same for all of the jets in our sample, and 
the differences in the physical conditions of the gas are due 
to the fact that the HH 111 and HH 34 jets are denser. 

- The total density is quite low (n H ~ 10 3 cm -3 ). From these 
values we infer mass fluxes of ~ 8 10~ 8 M yr -1 . Note that 
these values are upper limits because we are not taking into 
account the filling factor of the emitting gas for these jets. 

- The linear momentum fluxes are of ~ 1-4 1CF 5 
M yr~' kms _I . 

The above results demonstrate the potential of our com- 
bined optical/NIR analysis to study the jet physical structure 
and represent a good basis to plan future higher spectral-spatial 
observations with VLT-CRIRES/UVES and AO or interfero- 
metric instruments. A higher spectral resolution, in fact, would 
allow us to analyse the kinematic structure of the jets, while 
with higher spatial resolution we could extend our analysis to 
the base of the jet, where it is launched and accelerated, testing 
the theoretical MHD models. 
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